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CN 1 Abstract 
O 

Data from Z decays in DELPHI have been searched for — > D* + £~z7^ with 
the D* + decaying to D°tt + and D° -> K~7r+, K~7r + 7r + 7r or KT7r + (7r°). These 
events are used to measure the CKM matrix element |V C £,| and the form factor 

V 1 slope, p\: 

Qh 

J> ■ ^d*(1) |V c6 | = 0.0392 ± 0.0018 ± 0.0023; p\ = 1.32 ± 0.15 ± 0.33 

corresponding to a branching fraction: 

BR(B° -> D*+ri7 £ ) = (5.90 ± 0.22 ± 0.50)%. 

Combining these and previous DELPHI measurements gives: 

^d*(1) | Veil = 0.0377 ±0.0011 ±0.0019, = 1.39 ± 0.10 ± 0.33 and 
BR(B° -> D* + £-J7 e ) = (5.39 ± 0.11 ± 0.34)%. 

Using T D *{1) = 0.91 ± 0.04, yields: 

iVcfrj = 0.0414 ± 0.0012(stat) ± 0.0021 (syst.) ± 0.0018(t/ieon/). 

The 6-quark semileptonic branching fraction into a D* + emitted from higher 
mass charmed excited states has also been measured to be: 

BR(b -> D*+X£~Vt) = (0.67 ± 0.08 ± 0.10)%. 
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1 Introduction 

The Cabibbo, Kobayashi and Maskawa (CKM) matrix element V c & is a parameter of 
the Standard Model and its value needs to be fixed by experiments. This parameter 
determines the decay rate of 6-hadrons because |V u f,| governs the other possible charged 
weak decay of 6-quarks and contributes only 1-2% to the total decay rate. The value of 
| V c t| cannot be measured directly and there are two decay processes for which theoretical 
uncertainties are expected to be under control [1]: the inclusive semileptonic decay of 
6-hadrons corresponding to the b — > c£~Vi transition and the exclusive decay channel 

— > D* + £~V£. The latter is used in the following analysis, the D* + 1 is reconstructed 
through its decay to D°7r + and the D° meson is isolated using three decay channels: 
K~7T + , K~it + it + it~ and KT7r + (7r°). 

This study benefits from the reprocessing of DELPHI data taken between 1992 and 
1995 through improved versions of the event reconstruction algorithms. As a consequence, 
the number of signal events has increased by more than a factor two over those reported 
in [2] using the same decay final states. An additional decay channel of the D° (— > 
K~7T + (7r )) has been analysed which provides another factor of two increase. There 
are also improvements on the D* mass and 6-meson energy reconstruction. The D* + 
signal is now narrower which gives a better isolation of the signal over the combinatorial 
background. As the main source of experimental systematic uncertainty originates from 
the contribution of D* + mesons emitted in the decay of excited charmed states produced 
in 6-hadron semileptonic decays, additional observables have been defined to control the 
level of this contamination in a better way. The evaluation of the remaining contamination 
from the double charm cascade decays (b — > D* + DX, D — > t~V^Y) benefits from recent 
measurements of their rates. 



2 Measurement of |V C &| from the decay B[} — > D* + £ 

The value of |V c b| is extracted by studying the decay partial width for the process 
B° — > D* + £~I7^ as a function of the recoil kinematics of the D* + meson [1]. The decay 
rate is parameterized as a function of the variable w, defined as the product of the 
four-velocities of the D* + and B° mesons. This variable is related to the square of the 
four- momentum transfer from the to the t~T>i system, q 2 , by: 

ml+ + m |o - <f 
w = (1) 

B d 

and its value ranges from 1.0, when the D* + is produced at rest in the B[j rest frame, to 
about 1.5. Using the Heavy Quark Effective Theory (HQET) [3], the differential partial 
width for this decay is given by: 



dw 487r 3 
where K{w) contains kinematic factors: 



<lr ai[VJ \(w)TUw), (2) 



/CM = mU™B ~ m D ,fV^l{w + l) 2 (l + 4 ™ 1 ^ (3) 

\ w + 1 (1 — r) z J 



1 Throughout this paper charge-conjugate states are implicitly included. 
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with r = m D */m B and JF d »(w) is a hadronic form factor. 

Although the shape of the form factor, J-d*{iv), is not known, its magnitude at zero 
recoil, corresponding to w — 1, can be estimated using HQET. It is convenient to express 
Td* (w) in terms of the axial form factor h Al (w) and of the reduced helicity form factors 
Hq and H±: 



D*[W 



h 



aAw 



H 2 



Hi + m 



1 _|_ 4-w 1— 2wr+r 2 
1 w+1 (1-r) 2 



(4) 



The reduced helicity form factors are themselves expressed in terms of the ratios between 
the other HQET form factors (h v (w) : h A2 (w), h A[i (w)) and h Al (w): 



with 



H±(w) 



Ri(w) 



H (w) = 1 + 
\/l — 2wr + r 2 



w 



1-r 



1-r 
IT 



[I-R2H 



w — 1 



w + 1 



Ri(w) 



and R.2(w) 



h A Jw) + rh A Aw) 



(5) 
(6) 

(7) 



h Al (w) " h Al (w) 

Values for Ri(w) and R2(w) have been measured by CLEO [4] using different models. 
The unknown function h Al (w) is approximated with an expansion around w — 1 [5] : 



h Al H = h Al (l) x [1 - Sp 2 Al z + (53p 2 Ai - 15)z 2 - (231p 2 Ai - 91)z 3 ] , 



(8) 



/w + l+V2' 



An alternative 



where p Ai is the slope parameter at zero recoil and z - 
parametrization, obtained earlier, can be found in [6]. 

In the heavy quark limit (m b — > 00), ^ 7 d*(1) — h Al (l) coincides with the Isgur-Wise 
function [7,8] which is normalized to unity at the point of zero recoil. Corrections to 
J-d*(1) have been calculated to take into account the effects of finite quark masses and 
QCD corrections They yield F D *{1) = 0.91 ± 0.04 [9]. 

Experiments determine the product ^ 7 d*(1) |V c &| by fitting this quantity and the slope 
using the expression (2), convoluted with the experimental resolution on the w 
variable. Since the phase space factor )C(w) tends to zero as w — > 1, the decay rate 
vanishes in this limit and the accuracy of the extrapolation relies on achieving a reasonably 
constant reconstruction efficiency in the region close to w — 1. 

Results of the following analysis are expressed in terms of the q 2 variable. 



3 The DELPHI detector 

The DELPHI detector and its performance have been described in detail elsewhere [10]. 

The tracking system consisted of the Vertex Detector (VD), the Inner Detector (ID), 
the Time Projection Chamber (TPC) and the Outer Detector (OD) in the barrel re- 
gion while the Forward Chambers (FCA, FCB) covered the end-cap regions. The average 
momentum resolution for high momentum charged particles traversing the DELPHI mag- 
netic field of 1.2 T was cr(p)/p 2 = 0.0006 (GeV/c) -1 in the polar angle region between 
30° and 150°. 

The VD surrounded the beam pipe and consisted of three concentric layers of silicon 
micro-strip detectors at radii 6, 9 and 11 cm. Until 1994 the VD layers were single-sided 
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and provided only information in the R<p plane 2 . In 1994 the innermost and outermost 
layers were replaced by double-sided silicon micro-strip modules providing both Rxp and 
Rz measurements. The ID was placed outside the VD and consisted of a jet chamber 
providing R(f) information, and a trigger chamber providing a measurement of the z 
coordinate. In 1995 a new ID was installed with a longer jet chamber and straw tubes 
replacing the trigger chambers. The VD and ID were surrounded by the TPC, the 
main DELPHI tracking device. The TPC provided up to 16 space points per particle 
trajectory between radii of 40 and 110 cm. The OD consisted of 5 layers of drift tubes 
and complemented the TPC by improving the momentum resolution of charged particles. 

Hadrons were identified using the specific ionization (dE/dx) in the TPC and the 
Cherenkov radiation in the barrel Ring Imaging CHerenkov detector (RICH) placed be- 
tween the TPC and OD detectors. The muon identification relied mainly on the muon 
chambers, a set of drift chambers giving three-dimensional information located at the pe- 
riphery of DELPHI after approximately 1 m of iron. Electron identification relied mainly 
on the electromagnetic calorimeter in the barrel region (High density Projection Cham- 
ber, HPC) which was a sampling device having a relative energy resolution of ±5.5% for 
electrons with 46 GeV/c momentum, and a spatial resolution of ±2 mm in z. 

4 Hadronic Event Selection and Simulation 

Hadronic Z decays collected by DELPHI between 1992 and 1995 have been analysed. 
Each event was divided into two hemispheres by a plane orthogonal to the thrust axis. 
To ensure that the event was well contained inside the fiducial volume of the detector, 
the cosine of the polar angle of the thrust axis of the event had to lie between -0.95 
and +0.95. Charged and neutral particles were clustered into jets by using the LUCLUS 
algorithm [11] with the resolution parameter dj i n = 5 GeV. 

About 3.4 million events were selected from the full LEP1 data sets. The JETSET 7.3 
Parton Shower program [11] was used to generate hadronic Z decays, which were followed 
through the detailed detector simulation DELSIM [10] and finally processed by the same 
analysis chain as the real data. A sample of about nine million Z — > qq events was used. 
To increase the statistical significance of the simulation, an additional sample of about 
3.6 million Z — > bb events was analysed, equivalent to about 17 million hadronic Z decays. 
Statistics for the hadronic samples are given in Table I. 



Year 


Real data 


Simulated 
Z^qq 


Simulated 
Z -> bb 


1992+1993 
1994+1995 


1355805 
2012921 


3916050 
5012881 


1096199 
2495335 


Total 


3368726 


8928931 


3591534 



Table 1: Analysed number of events. In 1992 and 1993 only two-dimensional vertex 
reconstruction was available. 



2 In the DELPHI coordinate system, z is along the electron beam direction, <j> and R are the azimuthal angle and radius 
in the xy plane, and is the polar angle with respect to the z axis. 
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5 Selection of events for analysis 

Events corresponding to candidates for the decay — > D* + l~Ve are selected by 
requiring the presence of an identified lepton and of a D* + candidate in the same event 
hemisphere which is defined by the direction of the jet containing the lepton. D* + are 
measured using the decays to D°7r + . D° mesons are reconstructed using their decays into 
K~7r + , K~7r + 7r + 7r~ and K~it + (-k ). 

To reduce the contribution from Z decays into light flavours, the standard DELPHI 
event 6-tagging variable [10] is used to enhance the sample in Z — > bb decays. The b- 
tagging variable is essentially the probability that the analysed event originates from 
light quarks and is required to be less than 0.5 for D° — > KTtx + and K _ 7r + (7r°) decays 
and less than 0.1 for the D° — > K~n + -n + n~ decay. In addition, the mass of the D* + — i~ 
system is restricted to the range between 2.5 and 5.5 GeV/c 2 . 

5.1 Lepton identification 

Muons and electrons with momentum larger than 2 GeV/c and at least one associated 
hit in the VD are selected. 

Muons are identified using standard algorithms [10] based on the matching of the track 
reconstructed in the tracking system to the track elements provided by the barrel and 
forward muon chambers. Loose selection criteria are applied and the efficiency is ~ 80% 
for ~ 1% probability of hadron misidentification. 

Electrons are identified using a neural network algorithm providing about 75% effi- 
ciency within the calorimeter acceptance. The probability for a hadron to fake an electron 
was about 1%. Electrons from photon conversions are mainly produced in the outer ID 
wall and in the inner TPC frame. About 80% of them were removed, with negligible loss 
of signal, by reconstructing the conversion vertex. 

5.2 Isolation of the D° — > K 7r + decay channel 

The kaon candidate corresponds to a particle with the same charge as the lepton, 
with a momentum larger than 1 GeV/c and not identified as a pion by the standard 
algorithms [10] which combine information provided by the ionization deposited in the 
gas volume of the TPC and by the RICH detectors. The pion candidate must have a 
charge opposite to the kaon and a momentum larger than 0.5 GeV/c. The pion and kaon 
candidates must both have at least one associated VD hit in R0 and be situated in the 
same event hemisphere as the jet containing the lepton. The two tracks must intersect 
in space to form a D° decay vertex candidate and those with a x 2 probability lower 
than 10 -3 are rejected. K~7r + systems with a mass between 1.81 and 1.92 GeV/c 2 and 
with a momentum larger than 6 GeV/c are selected as signal candidates. Resolutions on 
the reconstructed D° mass measured on real and simulated events are given in Table 2. 
The selected mass window for the signal corresponds to about ±3(7. A D° track is then 
reconstructed using the parameters of the K~ and n + tracks fitted at their common vertex 
and imposing the condition that the D° mass, quoted as m(D°) in the following, is 1.8645 
GeV/c 2 [9]. 

The B decay vertex is obtained from the intersection of the D° and the lepton trajec- 
tories. This vertex must have a x 2 probability larger than 10~ 3 and only D° — 1~ pairs of 
total momentum larger than 10 GeV/c are kept. The B decay vertex is then required to 
be at a minimum distance from the position of the beam interaction point. This algebraic 
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D u decay channel 


92-93 MC 
MeV/c 2 


92-93 data 
MeV/c 2 


94-95 MC 
MeV/c 2 


94-95 data 
MeV/c 2 




16.0 ±0.8 


19.8 ±2.0 


16.7 ±0.6 


20.3 ± 1.4 


K _ 7T + 7r + 7T~ 


10.4 ±0.7 


12.7 ± 1.8 


10.5 ±0.5 


12.1 ± 1.0 



Table 2: Mass resolution for D° signals measured on real and simulated events. 



distance is evaluated along the direction of the D° — £~ momentum. A minimum distance 
between the D and B measured decay points, evaluated along the D° momentum, is also 
required. These conditions depend on the number of z VD hits associated to the tracks. 
If the decay distance is not measured along z with the VD, only cuts on decay distances 
transverse to the beam direction are applied. These requirements, which are rather loose, 
are given in Table 3. 



D u decay channel 


B vert, -v^ main vert. 


D vert. B vert. 


K~7r + 7r + 7r~ 

K"7T + (7T ) 


£/a e > -1, -2 
Ijot > 1, 2 
11 at > 2, 2 


£/a e > -2, - 1 

e/<T t > -i, - 0.5 

£/a e >-l, -1 



Table 3: Minimum requirements on the decay distance (£) between the B decay and 
the main vertices and also between the D and the B decay vertices. The first number 
corresponds to the cut in space whereas the second is the transverse distance cut, which 
is applied only when the z coordinate is not measured. Negative distances corespond to 
positions behind the beam interaction point. 

The D* + signal is identified by its decay to D°7r + . Each particle of charge opposite to 
the lepton candidate and emitted in the same event hemisphere as the jet containing the 
lepton is considered as a candidate for the n + . The track of this particle must form a 
vertex with the D° and the charged lepton trajectories and the vertex fit probability has 
to be higher than 10~ 3 . Signals for the cascade decay D* + — > D°7r + correspond to a peak 
in the distribution of the mass difference 5m = m(D°7r + ) — m(D°). 

The global efficiencies to select signal events have been estimated using the simulation 
(see Table 4), accounting for all analysis steps described above, apart from the branching 
fractions of the D* + and of the D° into the selected decay channels. 



D u decay channel 


92-93 MC 


94-95 MC 


K-7T+ 


(19.2 ±0.7)% 


(22.3 ±0.6)% 


K - 7T + 7r + 7r~ 


(8.6 ±0.3)% 


(10.8 ±0.3)% 


K-7T+(7T U ) 


(8.7 ±0.3)% 


(10.4 ±0.2)% 



Table 4: Global efficiencies of the analysis chain to reconstruct and select simulated 
signal events. The quoted uncertainties are statistical. 
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5.3 Isolation of the D° — > K 7r + 7r + 7r decay channel 

Similar selection criteria to those which were applied to isolate the D° — > KTit + decay 
channel, are used. Differences in the algorithm are related to the final state multiplicity. 
Each of the three pion candidates must have a momentum larger than 0.5 GeV/c and 
the total charge of the three pion system has to be opposite to the kaon charge. At least 
two, among the four charged particle track candidates for D° decay products must be 
associated to at least one VD hit in R0 and be situated in the same event hemisphere 
as the jet containing the lepton. As the reconstructed signal is narrower, KT-k + -k + -k~ 
systems with a mass between 1.84 and 1.90 GeV/c 2 are selected. The measured mass 
resolutions are given in Table 2. Because of the higher combinatorial background, as 
compared with the D° — > K _ 7r + decay channel, cuts on algebraic distances between the 
B and the main vertex and between the D and the B decay vertex are more severe than 
for the previous channel. They are summarized in Table 3. 

The same set of four particles can give two KT7r + 7r + 7r~ mass combinations if there is 
an ambiguity in the definition of the K~ and 7T~ candidates. Only one combination is 
kept in the analysis by using criteria which are based on the available particle identifica- 
tion information provided by the RICH and the TPC or, if this information is missing, 
assuming that the K~ has the larger momentum. 

The same selection criteria, as for the decay D° — > KT7r + , are applied to search for a 
D* + signal. The global efficiencies to select signal events have been estimated using the 
simulation (see Table 4), accounting for all analysis steps described above, apart from 
the branching fractions of the D* + and of the D° into the selected decay channels. 

5.4 Isolation of the D° — > K 7r + (7r°) decay channel 

The same criteria are applied, as in Section 5.2, to select the K~ and tt + candidates 
apart from the cut on the K~n + mass which is required now to be between 1.5 and 
1.7 GeV/c 2 . This mass interval corresponds to the satellite peak position for the decay 
D° — > K~p + when the 7r° emitted from the p + is soft. An estimate of the 7r° 4-vector is 
obtained by assuming that the decay is of the type D° — > K~p + , p + — > 7r + 7r° and the D°, 
p + and 7T° masses are used as constraints. In addition it has been assumed that the ir° is 
contained in the plane defined by the K~ and the ir + . When two solutions are possible, 
a choice is made according to criteria which have been defined using simulated events. 

The D* + signal is identified by its decay to D°7r + . Each particle of charge opposite 
to the lepton candidate and emitted in the same event hemisphere as the jet containing 
the lepton is considered as a candidate for the 7r + . The track of this particle must form 
a vertex with the D° and the charged lepton trajectories and the vertex fit probability 
has to be higher than 10~ 3 . Signals for the cascade decay D* + — > D°7r + correspond to a 
peak in the distribution of the mass difference 5m = m(D°7r + ) — m(D°). The peak is 
broader than for cases in which the D° was completely reconstructed using its charged 
decay products. Cuts on decay distances between the primary, the B and the D vertex 
are given in Table 3. 

The global efficiencies to select signal events have been estimated using the simulation 
(see Table 4), accounting for all analysis steps described above, apart from the branching 
fractions of the D* + and of the D° into the selected decay channels. The event selection 
described above does not ensure that only D° decaying into the K~7r + 7r° channel are 
selected. The simulation predicts that about 67% are of this origin and that there are also: 
K~£ + i/£ (18%), K~7r + X (3%) where X corresponds to neutrals, K _ K + (3%) where the 
K + is assumed to be a ir + and the remaining 10% originates from various other channels. 
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Apart from the last contribution, efficiencies have been determined for each individual 
channel using the simulation. Table 4 shows the selection efficiency corresponding to the 
weighted average for these channels. The branching fractions measured for each of these 
channels have been used for the real data (apart from K~7r + X which is assumed to be the 
same as in the simulation and equal to 5.6%, with an error of 0.6%) and a corresponding 
effective efficiency has been evaluated. A correction factor for the remaining 10% of 
the events of undetermined origin has been obtained from a fit to the simulation using 
the efficiencies of the four identified contributions (see Table 4) and ensuring that the 
simulated BR(B[j — > D* + £~Ve) is recovered. This correction has been used on real data 
with a relative uncertainty of ±25%, corresponding to the statistical error of the fit on 
simulated events. 

5.5 Selected event candidates 

The mass difference distributions corresponding to the variable 5m = m(D°7r + ) — 
m(D°) obtained for the three channels, are shown in Figure 1. 

The numbers of D* candidates obtained by fitting these distributions with a Gaussian 
(D° — ■» K~7r + and K _ 7r + 7r + 7r~) or a gamma distribution (D° — > K~7r + (7r )) for the signal, 
and a smooth distribution for the combinatorial background 3 are given in Table 5. 



Data set 


92-93 


94-95 


K _ 7T + 

KT7r + 7r + 7r~ 

K-7T+(7T ) 


193 ± 15 
144 ± 14 
286 ± 24 


328 ± 16 
243 ± 17 
494 ± 27 



Table 5: Number ofD*—£ candidate events selected in the two data taking periods and 
for the three D° decay channels. 



5.6 q 2 measurement 

As explained in Section 2, to measure |V C &| it is necessary to study the q 2 dependence 
of the differential semileptonic decay partial width dr(B[J — > D* + £~Vi)/dq 2 . For signal 
events, corresponding to the semileptonic decay B[J — > D* + £~V e , the value of q 2 has been 
obtained from the measurements of the and D* + four-momenta: 

q 2 = {Pi + Vvtf = (p^o-Pd*+) 2 - (9) 

d 

The D* + 4-momentum is accurately measured, as all decay products correspond to re- 
constructed charged particle trajectories. 4 To improve on the determination of the B° 
momentum, information from all measured 6-decay products is used, including an eva- 
luation of the missing momentum in the jet containing the lepton and the positions of 
the primary and of the secondary vertex, which are used as constraints to define the di- 
rection of the 6-hadron momentum. The nominal B[j mass is also used as a constraint in 
this fit. The missing momentum in each jet has been evaluated by comparing the recon- 
structed jet momentum with the expectation obtained by imposing energy-momentum 

3 The distribution selected for the combinatorial background is bs m (6m) = (6m — m 7! ) a ° (X)fc=i a fe<5 mfc 1 )i with n = 
2 or 3 and ao = 0.5. 

4 For the D° — > K~7r+(7r°) decay channel the accuracy is reduced by about 10% because of the missing n°. 
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conservation on the whole event. Finally, a momentum dependent correction is applied 
to the reconstructed 6-hadron momentum so that it remains, for simulated signal events, 
centred on the generated value. 

The smearing of the q 2 variable is studied with simulated signal events. The function 
7Z(q 2 — q 2 ,q 2 ) gives the distribution of the difference between the values of the recon- 
structed q 2 , q 2 , for events generated with a given value q 2 . Twenty slices in q 2 s of the 
same width have been considered. Within each slice, TZ(q 2 — q 2 , q 2 ) is parametrized as 
the sum of two Gaussian distributions (see Figure 2). The two central positions of the 
Gaussians, their standard deviations and the fraction of events corresponding to the nar- 
rower Gaussian are parametrized with a linear dependence on q 2 s . Such parametrizations 
are obtained independently for two sets of ten slices. Typical values of these parametriza- 
tions correspond to q 2 resolutions of 0.3 and 2 GeV 2 with about 50% of the events included 
in the narrower Gaussian. Resolution distributions obtained for D° reconstructed with 
only charged particles and for the KT7r + (7r°) decay channel are compared in Figure 3. 

The cuts applied to select the events which require a minimum momentum on the 
lepton, the D* + and the D* + — i system and the cut on the minimum value for the 
mass of the D* + — i system can possibly introduce a bias in the q 2 s distribution. A q 2 
dependent acceptance correction, e(q 2 ) has been evaluated by comparing the simulated q 2 
distributions for signal events before and after applying all analysis cuts. This correction 
has been normalized such that it does not change the number of accepted events for 
which an overall efficiency has already been determined. The corresponding distribution 
is given in Figure 4. It is uniform and does not show evidence for any significant bias. A 
linear dependence for the acceptance gives: 

e( q 2 ) = (0.985 ± 0.026) + (0.0024 ± 0.0043) x q 2 , q 2 in GeV 2 , (10) 

which is compatible with unity within quoted uncertainties. 

As the cuts used in the analysis are very similar for all data samples, the same q 2 
dependent acceptance correction has been used for all channels and data samples. 

6 The analysis procedure 

The purpose of this analysis is to determine the values of the parameters JF D *(1) | V cb | 
and p 2 Ai introduced in Section 2, using the measured q 2 distribution of candidate events. 
The predicted q 2 distribution for the signal is obtained using the theoretical distribution 
corrected by the overall efficiency and the q 2 s dependent acceptance, and then convoluted 
with the expected resolution function TZ(q 2 — g 2 ,? 2 )- The q 2 distributions for the other 
event sources are taken from the simulation or from the real data for the combinato- 
rial background. The q 2 distributions are rather similar for the signal and other event 
categories because the procedure used to evaluate q 2 from the B° and D* + 4-momenta 
overestimates the real q 2 value for background events. This is because algorithms have 
been defined for signal events and thus do not include the additional hadrons emitted 
in background sources. To enhance the separation between the signal and other event 
sources, three other variables have been used. As a result, the branching fraction for D* + 
production in the decay of higher mass charm states, BR(6 — > D* + Xl~Ve), has also been 
measured. Before describing these quantities, the different event classes contributing to 
the analysis are explained. 
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6.1 The Event Sample composition 

In addition to the signal (Si), which corresponds to the decay B° — > D* + £~V£, there 
are six classes of events which contribute to the background: 

• the combinatorial background (B) under the D* + peak; 

• real D* + — £~ events with the D* + produced in the decay of an excited charmed state 
(S2). These events correspond to the decay chain b — > D**£~V£, D** — > D* + X. In 
the present analysis, D** includes resonant as well as nonresonant Dnn systems; 

• real D* + — t~ events with the lepton originating from the decay of another charmed 
hadron (S3); 

• events in which the D* + is emitted during the hadronization of a charmed quark jet 
in Z — > cc events (S 4 ); 

• Z — > bb events with a real D* + candidate accompanied by a fake lepton of opposite 
sign (S 5 ); 

• real D* + — t~ events with the lepton originating from the decay of a r~ lepton (S 6 ). 

6.2 Separation of signal from background events 

There are two main classes of events which either do or do not contain a real D* + . The 
variable 5m (= m(D°7r + ) — m(D )) allows the two classes to be separated (see Figure 1). 
Variables, d±, are used to separate the different classes of events with a real D* + . They 
are obtained from a measurement of the number of charged particle tracks (excluding 
the charged lepton, the pion coming from the D* + and the D° decay products) which 
are compatible with the 6-decay vertex or with the main vertex. For the signal (Si), 
it is expected that all other charged particles in the 6-jet are emitted from the beam 
interaction region. This will be also true for (S4) , the remaining background from Z — > cc 
events, and for (Sq). For the other classes (S2, S3 and S5) it is expected that, for most of 
the events, one or more additional charged particles are produced at the 6-vertex. 

The variables d± are defined in the following way: 

• all charged particles, other than the D* + decay products and the lepton, emitted in 
the same event hemisphere as the 6-candidate, with a momentum larger than 500 
MeV/c, which form a mass with the D* + — i~ system lower than 6 GeV/c 2 and 
which have values for their impact parameters to the 6-decay vertex smaller than 2 
and 1.5 a" in R0 and z respectively, are considered; 

• selected particles, having the same (+) or the opposite (— ) charge as the lepton are 
considered separately. If there are several candidates in a class, the one with the 
largest impact parameter to the main vertex is retained and the quantity: 

x± = e(R0) x nsig(R(p) 2 + e(z) x nsig(z) 2 (11) 

is evaluated, where e and nsig are, respectively, the sign and the number of standard 
deviations for the track impact parameter relative to the main vertex. The sign 
of the impact parameter is taken to be positive (negative) if the corresponding 
track trajectory intercepts the line of the jet axis from the main vertex downstream 
(upstream) from that vertex. 

As the track impact parameters can extend to very large values because of the relatively 
long decay time of 6-hadrons, the variables d± are taken to be equal to the logarithm of 
(1 + x±) and their sign is taken to be the same as x±. For events with no spectator track 
candidate, that is, with no additional tracks compatible with the 5-decay vertex, a fixed 
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value of -4. is used for d±. Examples of distributions of the variable d + for the signal 
and for the different background components, corresponding to all analysed channels, are 
given in Figure 5. 

Due to track reconstruction effects, only 77% of signal events have no spectator track 
candidate instead of the expected 100%. Similarly, for the D** background, 1/3 of events 
with no additional tracks are expected by isospin at the 6-vertex whereas 46% are ob- 
served. These values allow the probabilities V(0\0) and V(0\ ^ 0) to be extracted for 
getting no spectator candidate when, respectively, there is not and when there is really 
such a candidate at generation level. Values for these two probabilities are respectively 
equal to 77% and 30% with a spread of ±5% corresponding to the different years and 
channels. As it will be explained later (section 6.3.3), these two quantities have been used 
to correct the present simulation of double charm cascade decays as it includes only the 
channel B — ► D~DX whereas the other contributions (B — > D~DX and D°DX) have a 
different charged particle topology. 



6.3 Fitting procedures 

Six event samples have been analysed separately corresponding to different detec- 
tor configurations (1992-1993 and 1994-1995) and to different decay channels of the D° 
(K _ 7r + , K^n + n + 7i~ and KT7r + (7r )). The analysis procedure is explained for a single such 
sample in the following sections 6.3.1 to 6.3.7. It is applied to all samples simultaneously 
to obtain the measurements. 

For each event (i), four measurements have been used: Xi = (q%, 5m, d + , cL)j. 
The parameters Fd*{X) |V c (,|, p 2 Al and the background from D** decays (S2) are obtained 
by minimizing a negative log-likelihood distribution. Other parameters, given in the 
following, have to be introduced to account for the various fractions of contributing event 
classes and to describe their behaviour in terms of the variables analysed. The likelihood 
distribution is obtained from the product of the probabilities to observe X{ for each 
considered event. These probabilities can be expressed in terms of the corresponding 
probabilities for each event's class and of their respective contributions in the event 
samples analysed: 

B x b(xi) + Si x Sj(xi) 
Vm = 1 ; «~ — ■ (12) 

b + T5 =1 s s 

In this expression, B and Sj, (j = 1, 6) are the numbers of events (fitted) corresponding 
to the combinatorial background and to the different classes of events with a real D* + . 
The functions b(x) and Sj(x), (j = 1, 6) are the respective probability distributions of 
the variable x. Each probability distribution for the x variable is considered to be the 
product of four probability distributions corresponding to the four different variables. 

These distributions can be obtained from data (b(x)) or from the simulation. The 
fitting procedure involves minimizing the quantity: 

N evt 

-m£ = -^lnP(fi), (13) 

i=i 

where N evt is the total number of events analysed. 

From external measurements there are also constraints on the expected number of 
events corresponding to the categories S3-S6. These constraints can be applied assuming 
that the corresponding event numbers follow Poisson distributions with fixed average 
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values (Sj). This is obtained by adding to Equation (13) the quantity: 

6 6 

-5>>s° + J>[r(^ + i)]. (14) 

3=3 3=3 

A similar expression is also added to account for the fact that the total number of 
fitted events must be compatible with the number (N) of selected events: 

-N f \nN + \nr(N f + l), (15) 

in which Nf, the number of fitted events, is equal to: Nf = ^2 i=1 6 Si + B. 

The list of fitted parameters is given in the following for each component contributing 
to the event sample analysed. 

6.3.1 signal events 

s i,q%(Qr)'- this distribution results from the convolution of the theoretical expected dis- 
tribution dF( ' Bd ~^^ - — (corrected by the q 2 s dependent efficiency and ac- 
ceptance) with the resolution function 7Z(q 2 — q 2 ., q 2 ). It depends mainly on 
p\ and on the assumed q 2 dependence for the ratio Ri and R 2 between the 
different contributing form-factors. 
Si,<5m(^"7.) : is a Gaussian distribution corresponding to the D* + signal for the D° — > 
K~7r + or K~7r + 7r + 7r~ decay channels and a gamma distribution for D° — > 
K _ 7T + (7r°). The two parameters for each distribution have been obtained 
from a fit to data. 

s i,d±{d±)'- these distributions are obtained from simulated signal events. The two dis- 
tributions, for the d + and gL variables are rather similar with about 77% 
probability for having no spectator track candidate and the remaining 23% 

being concentrated around zero. 
Si: the number of signal events can be expressed as: 

Si = N H x R b x 4 x f w x BR(B° T>* + rv t ) x BR(D*+ -> D°tt + ) x BR(D° X) x e(X). (16) 

d 

In this expression, Nh is the number of hadronic events analysed (Table 1), 

Rb is the fraction of hadronic Z decays into bb pairs, the factor 4 corresponds 

to the two hemispheres and the fact that muons and electrons are used, f^o 

u d 

is the production fraction of B[J mesons in a 6-quark jet, BR(B[J — > D* + l~V t ) 
is the semileptonic branching fraction of B° mesons which is measured in this 
analysis, 5 the other two branching fractions correspond, respectively, to the 
selected D* + and D° decay channels (Table 10), and e(X) are the efficiencies, 
given in Table 4, of the cuts applied in the analysis to select signal events. 
Note that S\ is proportional to (JFd.(1) |V c b|) 2 . 

6.3.2 events from D** decays 

These are events from the S2 class corresponding to the cascade decay b — > 
D**£~V t , D** -> D*+X. 

5 It is the integral of Equation 2 (divided by the total B^j width) and depends on the two fitted quantities (1) |V C {,| 
and p 2 Ai . 
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s 2,q%(<lr) : this distribution is taken from the simulation. Its variation for different 
fractions of D** states has been studied (see Section 7.3.4 and Figure 9) and 
accounted for as a small systematic shift and error. 
S2,Sm{Sm): the same distribution is used as for signal, si^ m {8m). 

S2,d±{d±): as for the signal, these distributions are taken from the simulation. It has 
been verified that they are not dependent on the type of D** state which 
produced the D* + . There is a marked difference between s 2 ^ + and s 2> d-, the 
latter being rather similar to the corresponding distribution for signal events. 
S 2 : the number of expected events is fitted without imposing constraints from 
external measurements. 



6.3.3 double charm cascade decay lepton events 

These are events from the S3 class corresponding to the cascade decay b — > 
D*+DX, D -> £~V e Y 

s 3,q 2 (Qr)'- this distribution is taken from the simulation. 
S3,Sm(Sm): the same distribution is used as for signal, Sij m {8m). 

s 3,d±(d±): when there are spectator tracks, the distribution s 2 ,d + (d + ) (with d + > —4), 
is used. The expected fractions of events with no spectator tracks in the d + 
and oL distributions have been evaluated from the measured contributions 
of D°D*+, D~D*+ and DjD*+ events [12,13], with the D -> £~X branching 
fractions and topological decay rates for the hadronic states X, taken from [9]. 
For d + it is expected that (39 ±6)% of the events have no spectator track 
and for d_ this fraction is (39 ±4)%. These numbers have to be corrected for 
reconstruction effects using the variables P(0|0) and V(0\ 7^ 0) introduced in 
Section 6.2. 

S3: the expected number of events from this source is taken from present mea- 
surements of b — > DDX decay rates which correspond to: 

BR(b -> D* + tI) + BR(b -> D*-£ + X) = (0.83 ± 0.21)%, (17) 

where the lepton originates from the D — > Y semileptonic decay. This value 
has been obtained using measurements from ALEPH [12] and BaBar [13] 
on exclusive double charm decay branching fractions of 6-hadrons, with a 
charged D* emitted in the final state, and using the inclusive semileptonic 
decay branching fractions of charmed particles given in [9]. 

Simulated events contain double charm decays of the type b — > D* + D s X 
only, with a corresponding branching fraction: BR(6 — > D* + £~X) = 0.25%. 
This rate has been rescaled to correspond to the value given in Equation 
(17), assuming that the experimental acceptance is similar for the different 
contributing channels. 



6.3.4 Z — > cc events 

s 4,q'ti ( lr) : this distribution is taken from the simulation. 
SA,5m{.8m): the same distribution is used as for signal, Si^ m {8m). 
Si,d ± (d±): as for the signal, it is taken from the simulation. 

S4: the expected number of events from this source is taken from the simulation 
after having corrected for the small difference between the rates for D* + 
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production in c-jets between simulated and real events [14]: 

P(c -> D*+) = (0.2392 ± 0.0035) MC <-> (0.226 ± 0.0U) Data . (18) 

The remaining contamination from cc is expected to be very small (of the 
order of 1%). 

6.3.5 fake lepton events 

Only fake lepton events associated with a real D* + and not coming from cc events, 
have to be considered as the other contributions have been already included. 

s 5,g2(<?r) : this distribution is taken from the simulation. 
S5,s m {Sm): the same distribution is used as for signal, Si^ m (8m). 
s 5,d±(d±): as for the signal, it is taken from the simulation. 

S5: the expected number of events from this source is taken from the simulation 
after having applied corrections determined, using special event samples, to 
account for differences between the fake lepton rates in real and simulated 
data (see Section 7.3). 

6.3.6 semileptonic decays with a r 

These are events from the S 6 class corresponding to the cascade decay 

b -> d*+t-x, t- -> e-v e Y 

s 6,q 2 (Qr)'- ^is distribution is taken from the simulation. 

S6,6m($ m )'- the same distribution is used as for signal, Si^miSm). 

s 6,d±(d±): is the same as the signal distribution Si )< j ± (d±) 

the expected number of events from this source is obtained assuming that 
the production rate for 6-hadron r semileptonic decays is 0.223 ± 0.004 of the 
rate with a/iore [15]. As for the cc background, events from r decays are 
expected to give a small contribution, of the order of 1%. 

6.3.7 combinatorial background events 

Real data events are selected in the upper wing of the D* + mass peak between 
0.15 and 0.17 GeV/c 2 for D° — > K~n + or K~n + n + iv~ channels, and in the range 
0.17 - 0.22 GeV/c 2 for D° -> K-jv+in ). 

b q 2 (Qr)'- this distribution is taken from real data events located in the upper part of 
the Sm distribution. 

bs m (5m): the parametrization given in Section 5.5 has been used. The same mass 
dependence has been taken for the first four samples while a parametrization 
corresponding to different values for the coefficients has been obtained for 
D° — * K~7T + (7r°) events. Parameters of these distributions have been fitted, 
outside the global likelihood fit, to the Sm distributions corresponding to the 
events selected for the analysis. 
bd ± (d±): as for the g 2 distribution, the d± distributions for combinatorial background 
events are obtained from analysed events, selecting those situated in the 
upper part of the Sm distribution. 
B: in each of the six samples, the total number of combinatorial background 
events is fitted over the total Sm range. 
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7 Measurements of Td*^) |V c &| and p A 

The six event samples have been analysed in the same way. Efficiencies and probability 
distributions have been determined independently for each sample. Common parameters 
corresponding to the description of physics processes have been fitted or taken from 
external measurements. The central values and uncertainties used for the latter are 
summarized in Table 10. 

7.1 Results on simulated events 

Signal events generated using the DELPHI simulation program correspond to a given 
dynamical model, using a given modelling of the decay form factors. The generated q 2 
distribution has been fitted using a parametrization derived from the one given in Sec- 
tion 2. As the model used in the simulation is a priori different from HQET expectations, 
it has been necessary to add arbitrary terms in the expression so that the fit will be rea- 
sonable over the whole q 2 range. These terms correspond to a polynomial development 
in powers of (w — 1), starting with at least quadratic terms so that they have no effect 
on the slope nor on the absolute value of the spectrum at the end-point corresponding to 
w — 1. 

Using the total number of generated events to fix the normalization, the equivalent 
values for the two parameters defining the signal in the simulation: 

J~d* (1) |V cb | = 0.03552 ± 0.00016; p\ = 1.088 ± 0.021 (19) 

are obtained. The fitted semileptonic branching fraction is equal to: 

BR(B° -> D* + rv e ) = (5.091 ± 0.020)%, (20) 

which agrees with the exact value of 5.103% used to generate these events. 

The exercise is repeated on pure signal events using the reconstructed q 2 distribution. 
This predicted distribution now includes the effects of the experimental reconstruction of 
the q 2 variable and of the acceptance. This gives: 

•Fd'CL) |V c6 | = 0.03549 ± 0.00050; p\ = 1.119 ± 0.052. (21) 

The fitted semileptonic branching fraction is equal to: 

BR(B° -> D* + rv e ) = (5.004 ± 0.054)%. (22) 

Finally, using the sample of Z — > qq and bb simulated events, the signal parameters 
are determined, including the different background components giving (see Table 6): 

•Fd'(I) |V c6 | = 0.03579 ± 0.00063; p\ = 1.122 ± 0.061. (23) 

The fitted semileptonic branching fraction is equal to: 

BR(!J -> D* + rv e ) = (5.081 ± 0.065)%, (24) 

demonstrating that the fitting procedure gives the expected values for the sig- 
nal parameters correctly. The q 2 distribution for MC events selected within the 
5m = [0.144, 0.147] GeV/c 2 interval for the K~7r + and K~n + n + n~ channels and within 
the 5m = [0.14, 0.17] GeV/c 2 interval for K _ 7r + (7r°) is shown in Figure 6 with the 
contributions from the fitted components. 
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Data set 


•rD*{i-) \ V c b| 


Pa, 


nn/n0 . n*+^-T7 "\ 
DtX{D d — > U I Vl) {/o) 


ts. Ti yz-yo 
K-JT+ 94-95 


n 0^7^ -i- n nnon 
U.UorfO it u.uu/u 

0.0356 ± 0.0013 


1 97 -4- D 1 7 

1.16 ± 0.13 


r. i a _i_ n 91 

4.94 ± 0.14 


K^tt+tt+tt- 92-93 
K^tt+tt+tt- 94-95 


0.0356 ± 0.0020 
0.0363 ± 0.0014 


1.03 ± 0.21 
1.13 ± 0.13 


5.28 ± 0.23 
5.20 ± 0.15 


K"7r + (vr u ) 92-93 
K"7T+(7r°) 94-95 


0.0355 ± 0.0018 
0.0351 ± 0.0013 


1.14 ± 0.17 
1.05 ± 0.13 


4.95 ± 0.19 
5.06 ± 0.14 


Total sample 


0.03579 ± 0.00063 


1.122 ± 0.061 


5.081 ± 0.065 



Table 6: Fitted values of the parameters in Z — > qq and bb simulated events. The quoted 
uncertainties are statistical. 



7.2 Results from data 

To analyse real data events, additional corrections have been applied to account for 
remaining differences between real and simulated events. Central values and uncertainties 
on these corrections are explained in the following when evaluating systematic uncertain- 
ties attached to the present measurements. 

The results obtained on the six data samples and using the total statistics are given 
in Table 7. 



Data set 


^D*(l) |V cb 


Pi 


BR(B° -> D* + £-V e ) (%) 


K--K+ 92-93 
K~ti + 94-95 


0.0394 ± 0.0055 
0.0340 ± 0.0041 


1.15 ± 0.48 
0.71 ± 0.45 


6.55 ± 0.77 
6.11 ± 0.55 


K-n + n + n- 92-93 
K-tt+tt+tt- 94-95 


0.0410 ± 0.0058 
0.0342 ± 0.0042 


1.43 ± 0.46 
1.12 ± 0.41 


6.06 ± 0.77 
5.01 ± 0.51 


K-n+in") 92-93 
K-n + (7i ) 94-95 


0.0407 ± 0.0043 
0.0404 ± 0.0031 


1.36 ± 0.35 
1.48 ± 0.24 


6.22 ± 0.57 
5.70 ± 0.40 


Total sample 


0.0381 ±0.0018 


1.23 ±0.15 


5.83 ±0.22 



Table 7: Fitted values of the parameters on real data events. The quoted uncertainties 
are statistical. As 92-93 event samples have a reduced sensitivity to the D** background 
(S2), fitted values quoted in this Table, when corresponding to individual event samples, 
have been obtained using a fixed value for the branching fraction BR(6 — > D* + X£~ve) 
(=0.67 %, see Equation 29). Results quoted for the total statistics have been obtained 
letting free this quantity to vary in the fit. 

The values obtained are: 

^d*(1) |V C& | = 0.0381 ± 0.0018; p\ = 1.23 ± 0.15, (25) 
which correspond to a branching fraction equal to: 

BR(B° -> D* + rv e ) = (5.83 ± 0.22)%. (26) 

The correlation coefficient p(jFc,*(l) |V c t| , p 2 Ai ) is equal to 0.894. Fitted fractions of 
the different components are given in Table 8. 

Distributions of the q 2 and d± variables for events selected within the 
5m = [0.144, 0.147] GeV/c 2 interval for the K~ir + and KT7r + 7r + 7r~ channels and within 
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the 5m = [0.14, 0.17] GeV/c 2 interval for K 7r + (7r°) are shown in Figures 7 and 8 with 
the contributions from the fitted components. 



Signal 




Cascade 


Charm 


Fake lept. 




Comb. Backg. 


(Si) 


(S 2 ) 


(Sa) 


(S 4 ) 


(S 5 ) 


(s 6 ) T 


B 


1196 ±35 


319 ±38 


129 ± 11 


12 ±3 


67 ±8 


16 ±4 


523 ± 23 




26.7 ±3.2 


10.8 ±0.9 


1.0 ±0.3 


5.6 ±0.7 


1.3 ±0.3 


43.7 ± 1.9 



Table 8: Number of events and fitted fractions (in % of signal events) attributed to the 
different components of the analysed sample of events selected within the 5m mass interval 
corresponding to the D* + signal. 



7.3 Evaluation of systematic uncertainties 

Values for the parameters taken from external measurements and hypotheses used in 
the present analysis have been varied within their corresponding range of uncertainty. 
The results are summarized in Table 9. 



parameter 


central value 


rel. err. on 


rel. err. on 


rel. err. on 


or hypothesis 


and uncert. 


^.(l)|V c6 |(%) 


P% (%) 


BR(%) 


External parameters 










Rates and BR 


Table 10 and 11 


T2.5 


±0.4 


T5.3 


K~tt + X rates 


1.100 ±0.025 


T0.5 


0.0 


T-l.l 


b-hadron frag. 


see text 


T0.8 


T3.0 


0.0 


r(M) 


(1.542 ±0.016)ps 


T0.5 


0.0 


T-1.0 


Detector performance 










Tracking efficiency 


see section 7.3.2 


T 1.1 


0.0 


=F 2.2 


Lepton identification 


±1.5% (e), ±2.0% (/x) 


T0.5 


±0.1 


T-1.0 


Fake Lepton rates 


Table 12 


0.0 


±0.1 


0.0 


q 2 resolution 


see text 


±2.3 


±6.2 


T0.4 


q 2 acceptance 


see section 5 


±0.5 


±1.5 


T0.2 


Control of d± dist. 


see text 


±3.0 


±1.1 


±5.5 


Selection efficiency 


Table 4 


T0.6 


±0.1 


±1-2 


MC statistics 


see section 7.1 


±1.8 


±5.9 


±1.3 


Signal modelling 










Ri(w) and B.2(w) 


see text 


±1.0 


±22.8 


0.0 


Backg. modelling 










D** states 


see text 


±2.2 


±5.3 


±0.6 


Double charm cascade decay rate 


0.0083 ± 0.0021 


±0.4 


±1.7 


T0.4 


Bj D*+r-V T X 


0.0127 ±0.0021 


T0.4 


T0.4 


T0.6 


P(c-»D*+X) 


0.226 ±0.014 


0.0 


±0.1 


0.0 


Total systematics 




±5.8 


±25.2 


±8.4 



Table 9: Systematic uncertainties given as relative values expressed in %. The total 
systematics are obtained by summing the components in quadrature. 



7.3.1 Uncertainties related to external parameters 

• Values for D and D* branching fractions into the analysed final states, the R b value 
and the 6-hadron lifetime have been taken from [9]. A summary of the values used 
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parameter 
or hypothesis 


central value 
and uncert. 


BR(D*+ -> D°7T + ) 
BR(D° -> KTtt+) 

BR(D° -> K-TT + TT + TT-) 

BR(D° -> K-tt+tt ) 
BR(D° -> K-£+^) 
BR(D° -> K-X+) 


0.21664 ±0.00068 
0.388 ±0.013 
0.677 ±0.005 
0.0380 ±0.0009 
0.0746 ±0.0031 

0.131 ±0.009 
0.0686 ±0.0030 
0.00412 ±0.00014 



Table 10: Values for the external parameters used in the analysis. The quoted value for 
BR(D° — > K~t + V() corresponds to the sum of the branching fractions for the electron and 
muon final states. 



in the present analysis is given in Table 10. They have been varied within the corre- 
sponding range of uncertainty and the systematic errors induced in JF D »(1) |V c b] (%), 
p 2 Al (%) and the BR(B° -> B* + £~u e ) are given in Table 11. 

Global efficiencies to select D° — > KTTr+(Tr ) events have been estimated using the 
simulation as described in section 5.4. Measured branching fractions of several 
D° — * K~n + X decay channels have been used for real data and a correction factor 
has been applied to take into account events from undetermined origin. 
Simulated events have been generated using the JETSET 7.3 program with the par- 
ton shower option [11]. The non-perturbative part of the fragmentation of 6-quark 
jets is taken to be a Peterson distribution which depends on a single parameter, e^: 

N 

D(z) = 2. (27) 

L 2 1 — zJ 

In this expression, N is a normalization factor and z = -j^jp^f with B and 
b indicating the B hadron and the 6-quark. The average fraction of the beam 
energy taken by weakly decaying 6-hadrons has been evaluated in [14] to be 
< X E >= = 0.702 ± 0.008. Simulated events, generated with the param- 

eter = 0.002326, correspond to < X E >= 0.7035. The effects of a variation of 
the average value < X E > and of the shape of the fragmentation distribution, on 
the results of the analysis have been studied by weighting events generated with a 
known value of the variable z so that they correspond to recent measurements [16]. 
In addition to a change in the slope of the B momentum distribution, < X E > 
increases by 2%. 

A new parametrization of the resolution function TZ(q 2 — q 2 , q 2 ) described in section 
5.6 has been determined on simulated weighted events and the analysis has been 
redone giving relative variations of —0.8% and —3%, respectively, on JF D »(1) |V C &| 
and p 2 Ai . 

The 6-hadron lifetime used in the simulation is equal to 1.6 ps and is independent 
of the type of produced 6-hadrons. Events have been weighted so that their lifetime 
agrees with present measurements. Efficiencies given in Table 4 have been deter- 
mined using weighted events and uncertainties related to the present accuracy on 
the B° lifetime measurement can be neglected. 
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parameter 


rel. err. on 


rel. err. on 


rel. err. on 


01 nypomesis 


r n n ^ iv j 




Dixy /0 J 


R, 






-(-U.OZ, 


P(6-Bj) 


Tl.69 


±0.02 


^3.39 


BR(D*+ -> D°7T + ) 


^0.37 


±0.03 


^0.76 


BR(D° -> K"7T+) 


T0.31 


±0.01 


^0.63 


BR(D° -> K^TT + TT + TT-) 


^0.68 


0.0 




BR(D° -> K-7T+(7r )) 


Tl.70 


±0.35 


T3.64 


Total 


^2.5 


±0.4 


^5.3 



Table 11: Systematic uncertainties from external parameters as relative values expressed 
in %. The uncertainty related to the BR(D° — > iT~7r + (7r )) includes the contribution of 
the D° — > K^n + n°, D° — > K^i^vg and D° — > K~K + branching fractions. Central values 
and errors of these parameters are given in Table 10. 

7.3.2 Uncertainties from the detector performance 

• Differences between simulated and real data events on the tracking efficiency have 
been studied in [17] and correspond to ±0.3% for each charged particle. For the 
soft pion coming from the D* + an uncertainty of ±1% has been assumed. 

• Differences between simulated and real data events on lepton identification have 
been measured using dedicated samples of real data events [18] and the real data to 
simulation ratios are equal to (88.5±1.5) 92 _93% and (94.0±1.5) 9 4_ 95 % for electrons. 
For muons, the real data to simulation ratios are equal to 96% for the two periods 
with a ±2% uncertainty. 

• Differences between fake lepton rates in which the lepton is a misidentified hadron, 
have also been measured using dedicated samples of real data events and compared 
with the simulation [18] to obtain correction factors which are summarized in Ta- 
ble 12. 



Data set 


electron 


muon 


92-93 
94-95 


0.69 ±0.03 
0.77 ±0.03 


1.44 ±0.03 
1.61 ±0.03 



Table 12: Correction factors to apply to simulated events in which the candidate lepton 
is a misidentified hadron. 

• Resolution of the q 2 variable. 
A resolution function, common to all three D° decay channels has been used. This 
function is determined independently for the 92-93 and 94-95 data samples. To 
quantify the importance of controlling the experimental resolution on q 2 the widths 
of the fitted Gaussians have been increased by 5%. This value is two times larger 
than observed differences between the averaged missing energy measured in jets 
for real and simulated events. It corresponds to the increase in smearing of the 
resolution function when including events with a missing 7r°. The induced variations 
in JF D *(1) \V c b\ and p 2 Ai are ±0.3% and ±1.2% respectively. 

The uncertainty on the parametrization of the resolution distributions has been 
evaluated by varying the number of fitted groups of slices in q 2 on which a linear 
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variation of the parameters of the two Gaussian distributions were evaluated. Re- 
sults obtained with two groups of 10 slices and with five groups of four slices have 
been compared. This corresponds to relative variations on JF D , (1) |V C &| and p 2 Al of 
±2% and ±5% respectively. 

Results obtained when including or excluding K _ 7r + (7T°) events, which have a poorer 
resolution, in the determination of the resolution function have been also compared. 
This corresponds to relative variations on JF D ,(1) |V c6 | and p 2 Ai of ±1.0% and ±3.5% 
respectively. 

Measured differences obtained from these comparisons have been summed in 
quadrature. 

The value of q 2 is obtained from the measurements of the B and D* + 4-momenta 
(see Section 5.6). The B momentum is obtained from a constrained fit, imposing 
the B meson mass, which includes information from primary and secondary vertex 
positions and from the energy and momentum of the particles belonging to the jet 
that provide an estimate of the B momentum and direction. Uncertainties on the 
polar and azimuthal angles giving the B direction, and on the magnitude of the 
B momentum, which were determined from the measurement of the D* + , charged 
lepton and missing-jet momenta, have been varied by ±30% and new resolution 
distributions for q 2 have been obtained. Corresponding variations on fitted values 
for J^D'iX) |V c f,| and p Ai are found to be negligible. 

• Control of the d± distributions. Distributions of the d± variables obtained for 
events selected for values of 5m higher than the D* + signal, in real and simulated 
events, have been compared (see the bottom two distributions in Figure 5). The 
probabilities for having no spectator track differ by (2.5 ± 1.0)% between data and 
the simulation. To account for this difference the corresponding probabilities for no 
spectator track have been varied by ±3%, simultaneously for signal and background 
components with a real D* + . Such a variation does not apply for events from the 
combinatorial background as the shape of the corresponding distributions has been 
taken from real events. 

The effect of a different shape of the d± distributions has also been evaluated for 
double charm cascade decays (S3). A flat distribution has been considered for 
d± > 2 to account for the different topologies of D D* + , D~D* + and D~D* + . The 
effect of this variation has been found to be negligible. 

• The effect of a possible difference between the tuning of the 6-tagging [10] between 
real and simulated data events has been neglected because loose criteria have been 
used in this analysis. 

7.3.3 Uncertainties on signal modelling 

These uncertainties correspond to the use of the w dependent ratios R\(w) and R2{w) 
defined in Equation (7). Values for these quantities, using different models, have been 
obtained by the CLEO collaboration [4]: R l = 1.18±0.30±0.12, R 2 = 0.71±0.22±0.07 
with a correlation p(Ri,R 2 ) = —0.82 between the uncertainties on these two mea- 
surements. Relative variations induced in the fitted parameters JFd*(1) |V C &|, p 2 Al and 
BR(B° — > D* + £~Ve) have been obtained varying the values of R\ and R 2 within their 
corresponding range of uncertainty. They are given in Table 13. 

As observed already in previous analyses, the uncertainty on R 2 dominates the sys- 
tematic uncertainty on p 2 Ai . 
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AR 1 
AR 2 



±la 

±1(7 



AF D *(l)\V cb \ 
^D*(l)|V c6 | 



(%) 



T 1.3 
T 1-8 



± 0.3 
=F 22.5 



abr(b°->d*+^-i^) 

BR(Bf-^D*+^-17 £ ) 



(%) 



0.0 
0.0 



Table 13: Relative variations of the fitted parameters due to the R\ and R 2 measurements. 
The la variation corresponds to the sum in quadrature of the statistical and systematic 
uncertainties. 



7.3.4 Uncertainties on background modelling 

• The fraction of D* + mesons originating from decays of D** mesons depends on the 
total production rate of these states and on their relative fractions. 
Combining present measurements, the production rate of D* + mesons originating 
from D** decays and accompanied by an opposite sign lepton is [14]: 

BR(6 -> D* + Xrv e ) = (0.8 ± 0.1)%. (28) 

This information is not included in the fit as D* + events produced in D** decays 
are directly fitted, simultaneously with ^73* (1) |V C &| and p\ v to the data giving: 

BR(6 -> D* + Xrv £ ) = (0.67 ± 0.08 ± 0.10)%, (29) 

which is compatible with the expectation given in Equation (28). The statistical 
error correlation coefficients are: p(J r o*(l) |V c b| , BR(6 — > T)* + X£~T>g)) = —0.171 
and p(p 2 Ai ,BR(b — > D* + X£~Ve)) = —0.061. The quoted systematic, in Equation 
(29), has been evaluated by considering the same sources of errors as are listed in 
Table 9. 

To evaluate the effect of the uncertainty in the sample composition of produced D** 
states, the model of [19] has been used. Parameters entering into this model have 
been varied so that the corresponding production rates of the narrow states remain 
within the ±1<t measured ranges defined in Equations (30, 31): 



(0.63 ±0.10)% (30) 
(0.23 ± 0.08)% or < 0.4% at the 95% CL (31) 

0.37 ± 0.14 or < 0.6 at the 95% CL, (32) 

where R** is the ratio between the production rates of DJ; and D x in b- meson semilep- 
tonic decays. 

A dedicated simulation program has been written to generate the decay distributions 
of the different D** states. Correlations between the lepton and hadron momenta 
induced by the decay dynamics are included. The w dependence of the different 
form factors has been parametrized according to the model given in [19]. It has 
been assumed that, in addition to narrow states whose production fractions are 
given in Equations (30-32), broad D*7r final states, emitted in a relative S wave, are 
produced. 

The two sets of model parameters giving the two most displaced central values for 
the q 2 distribution are used to evaluate the systematic uncertainty coming from the 



BR(B — > T>i£~i) = 
BR(B — > Y)* 2 l~Vi) = 
BR(B -> D* 2 £~i7i) 
BR(B -»• D^-Vi) 
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sample composition of D** states . These two distributions are shown in Figure 9 
and the fitted values obtained with these two models are given in Table 14. 





Model 1 


Model 2 


•Fd*(1)|V c6 |(%) 
Pa, 

br(b° -> n*n-i7i)(%) 


0.0400 ±0.0017 
1.38 ±0.13 
5.92 ±0.21 


0.0383 ±0.0018 
1.25 ±0.15 
5.84 ±0.22 



Table 14: Fitted values corresponding to the two models describing D** production. 

The average of these two results is used to determine the central values for 
JF D »(1) |V C (,| and p 2 Al and half of their difference is taken as systematic uncertainty; 
this gives: 

Fd*(1) \Vcb\ = 0.0392 ± 0.0018; p\ = 1.32 ± 0.15. 

• The rate for the double charm cascade decay background, evaluated from simulated 
events, has been rescaled to agree with present measurements (see Equation (17)). 

• The small components of tau and charm backgrounds have been evaluated using 
present measurements. 

• The modelling uncertainty of the combinatorial component corresponding to events 
situated under the D* + peak has a negligible contribution. 

8 Combined result 

The present measurement of Fd*(X) \Vd>\ an d Pa, nas been combined with the previous 
DELPHI result [20] obtained with a more inclusive analysis in which the D* + was recon- 
structed using the charged pion and tracks attached to a secondary vertex, accounting 
for the D° decay products. The values obtained in the previous DELPHI analysis were: 

.F D .(1) |V c6 | = 0.0355 ± 0.0014^^3. p ^ = L34 ± 0.14 + °;^ and 
BR(B° -> D*+e-u t ) = (4.70 ± 0.13 

Modifying the central values and uncertainties of the parameters entering in this analysis 
so that they correspond to the values taken for the present measurement yields the results: 

^D'(l) |V c6 | = 0.0372 ± 0.0014 ± 0.0025; p\ = 1.51 ± 0.14 ± 0.37 and 
BR(B° -> D*+£-ui) = (5.20 ± 0.14 ± 0.42)%. 

The average with the present measurement has been obtained using the method 
adopted by the LEP Vcb working group [14]. Common sources of systematic uncer- 
tainties between the two analyses have also been identified and properly treated in the 
averaging procedure. The statistical correlation between the two measurements has been 
evaluated to be 8% (it was evaluated to be 5% in [2] where a similar combination was 
done). 

Averaging the two measurements gives the following results: 

^b'(l) |V c6 | = 0.0377 ±0.0011 ±0.0019; p% = 1.39 ± 0.10 ± 0.33 and 
BR(B~° -» D* + £-V e ) = (5.39 ± 0.11 ± 0.34)%. 

6 Values of the parameters (sec [19]) corresponding to Model 1 are: r' = —0.2, r(l) = 0.5, fi = —0.375 and T2 = 0.375. 
The corresponding values for Model 2 are: t' = —2.0, r(l) = 0.83, fi = 0. and T2 = 0.75. 
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Using F D *{1) = 0.91 ± 0.04 gives: 

|V cb | = 0.0414 ± 0.0012 ± 0.0021 ± 0.0018, 
where the last uncertainty corresponds to the systematic error from theory. 

9 Conclusions 

Measurements of JF D *(1) \V c b\, p 2 Ai and of BR(B[J — > D* + £~z7^) have been obtained 
using exclusively reconstructed D* + decays by the DELPHI Collaboration. Variables 
have been defined which allow different decay mechanisms producing D* + mesons in the 
final state to be separated. 

The following values have been obtained: 

•Fd* (1) |V c6 | = 0.0392 ± 0.0018 ± 0.0023; p\ = 1.32 ± 0.15 ± 0.33, 
which correspond to a branching fraction: 

BR(B° -> D*+£-V e ) = (5.90 ± 0.22 ± 0.50)%. 

These values are in agreement with previous measurements obtained by the 
ALEPH [21], DELPHI [20] and OPAL [22] collaborations. 

The 6-quark semileptonic branching fraction into a D* + emitted from higher mass 
charmed excited states has also been measured to be: 

BR(b -> D*+X£-V e ) = (0.67 ± 0.08 ± 0.10)%. 

Combining the present measurements with the previous analysis from DELPHI [20] 
which was done using a more inclusive approach, yields: 

J~d* (1) |V c6 | = 0.0377 ± 0.0011 ± 0.0019; p% = 1.39 ± 0.10 ± 0.33 and 
BR(B° -> D*+e-u e ) = (5.39 ± 0.11 ± 0.34)%. 

Using JFd*(1) = 0.91 ± 0.04 this gives: 

|V cb | = 0.0414 ± 0.0012 ± 0.0021 ± 0.0018, 

where the last uncertainty corresponds to the systematic error from theory. 
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Figure 1: 5m = m(D°7r + ) — m(D°) distributions for the D° — ► if _ 7T + (upper), 
D° — > if~7r + 7r + 7r~ (middle) and D° — > if~7r + (7r ) (lower) decay channels. Combina- 
tions with the wrong K-lepton charge correlation are superimposed as darker histograms. 
Events registered in 92-93 and 94-95 have not been distinguished. The curves show the 
fits to the right-sign distributions described in the text. 
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Figure 2: Fit of the TZ(q^ — q 2 r1 q^) slices for the 94-95 data taking period, as expected 
from simulated events. The simulated q$ central value is quoted above the corresponding 
slice. 
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Figure 3: Comparison between the resolution functions obtained for D° decay channels 
with and without a missing particle for the 92-93 and 94-95 data taking periods, as ex- 
pected from simulated events. 
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Figure 5: Distributions of the d + variable for signal and background components. All 
distributions have been normalized to unity. The content of the bin at d = —4. has been 
inserted on each plot, corresponding to events with no spectator track candidate. In the 
two lower plots, distributions obtained for combinatorial background events selected in 
real and simulated data can be compared. For real data the distribution for double charm 
cascade decays (S3) has been corrected as described in section 6.3.3. 
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Figure 6: Fit of MC qq and bb events. The three analysed D° decay channels and the two 
data taking periods have been combined. Only events selected within the Sm mass interval 
corresponding to the D* + signal are displayed. The small contribution of Z — >• cc decays 
and leptons originating from t~ events has been included in the fake lepton component. 
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Figure 7: Fit on real data events. All periods are combined. Only events selected 
within the 5m mass interval corresponding to the D* + signal are displayed. The small 
contribution of Z — ► cc decays and leptons originating from t~ events has been included 
in the fake lepton component. 
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Figure 8: Distributions for the d± variables, for events selected within the 5m interval 
of the D* + signal and corresponding contributions from the fitted components. The small 
contribution of Z — > cc decays and leptons originating from r~ events has been included 
in the fake lepton component. The measured (fitted) number of events in the first bin, 
at -4., are 1324 ± 36 (1358) and 1498 ± 39 (1540), for the d + and d_ distributions 
respectively. 
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Figure 9: q 2 distributions, normalized to unity, obtained using the two sets of parameters 
of the model [19], which correspond to the largest variation in the central value of these 
distributions and which give production rates for narrow D** states that are compatible 
with present measurements. The precise definition for models 1 and 2 is given in section 
7.3.4- The three components given in each histogram correspond, from top to bottom, to 
narrow 2^ , broad 1 + and narrow 1 + states. 



